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Are ring galaxies the ancestors of giant low surface brightness galaxies? 
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Giant low surface brightness galaxies (GLSBs), such as Malin 1, have unusually large and flat discs. Their formation is a 
puzzle for cosmological simulations in the cold dark matter scenario. We suggest that GLSBs might be the final product 
of the dynamical evolution of collisional ring galaxies. In fact, our simulations show that, approximately 0.5-1.5 Gyr after 
the collision which lead to the formation of a ring galaxy, the ring keeps expanding and fades, while the disc becomes 
very large (^ 100 kpc) and flat. At this stage, our simulated galaxies match many properties of GLSBs (surface brightness 
profile, morphology, HI spectrum and rotation curve). 
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1 Introduction 



The giant low surface brightness galaxies (GLSBs) are low 
surface brightness galaxies (LSBs) characterized by the un- 
usually large extension of the stellar and gaseous disc (up 
to ~ 100 kpc; Pickering et al. 1997 and references therein) 
and by the presence of a normal stellar bulge (Sprayberry et 
al. 1995; Pickering et al. 1997). Their prototype is Malin 1 
(Bothun et al. 1987). The existence of GLSBs is a puzzle for 
cosmology and in particular for cosmological simulations 
in the cold dark matter scenario. In fact, most cosmological 
simulations including a baryonic component produce discs 
which are too compact and too bulge-dominated to match 
the properties even of a Milky Way-like galaxy (D'Onghia 
et al. 2006). Thus, there is no way to form GLSBs, whose 
discs are flat and huge, in current cosmological simulations. 
Various mechanisms have been proposed for the origin of 
GLSBs, but none of them is able to completely solve the 
problem. For example, a large-scale bar can redistribute the 




Fig. 1 Projected density of stars of a simulated ring 
galaxy. Run A3 of Mapelli et al. (2008a), - 100 Myr af- 
disc matter and increase the disc scaldength(Noguchi2001).| the galaxy interaction. The frame measures 90 kpc per 

edge. The density, projected along the z— axis, scales lin- 
early from to 27 Mq pc~^. 



However, bar instabilities normally do not increase the disc 
scalelength by more than a factor of 2.5, which is not suffi- 
cient to produce the observed GLSBs. 

In this proceeding, we show that the propagation of the 
ring in an old collisional ring galaxy can lead to the redis- 
tribution of mass and angular momentum in both the stellar 
and gas component out to a distance of ~ 100 — 150 kpc 
from the centre of the galaxy, producing features (e.g. the 
surface brightness profile, the star formation, the HI emis- 
sion spectra and the rotation curve) which are typical of 
GLSBs. 



Corresponding author: e-mail: mapelli @physik.unizh.ch 



2 Methods 

We simulate galaxy interactions which lead to the formation 
of a collisional ring galaxy. The details about the simula- 
tions are reported in Mapelli et al. (2008a) and in Mapelli 
et al. (2008b). Here we remind that the simulations have 
been done with the N-body-SPH code GASOLINE (Wads- 
ley et al. 2004). Both the target and the intruder galaxy have 
a Navarro, Frenk & White (1996, NFW) dark matter halo. 
The target galaxy has a stellar and gaseous exponential disc 
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Mapelli & Moore: Ring galaxies and GLSBs 




Fig. 2 Projected density of gas and stars of a simulated 
ring galaxy. The density is projected along the z— axis. Run 
C of Mapelli et al. (2008b). Top left panel: initial conditions. 
Top right panel: ~ 100 Myr after the galaxy interaction. 
Bottom left panel: ~ 500 Myr after the galaxy interaction. 
Bottom right panel: ~ 1 Gyr after the galaxy interaction. 
The frames measure 260 kpc. The colour coding indicates 
the density, projected along the z-axis, in logarithmic scale 
(from 0.2 to TOM© pc-^). 



and a stellar bulge, and the ratio between the mass of the 
target and that of the intruder is ~ 2. 



3 Results 

The target galaxy develops a well-defined ring already ~ 
100 Myr after the interaction (Figure 1). The morphology 
of the simulated ring galaxy matches many of the proper- 
ties of the Cartwheel galaxy, including the 'spokes' (see 
the appendix of Mapelli et al. 2008a). The simulation also 
matches the surface brightness profile of stars in Cartwheel 
(Higdon 1995). 

We continue the simulation till ~ 1.5 Gyr after the in- 
teraction. We note that the ring-galaxy phase is quite short- 
lived: already ~ 0.5 Gyr after the interaction, the ring has 
propagated up to ~ 70 — 90 kpc and its surface density has 
significantly lowered (Fig. 2). At the same time, the stellar 
disc has became extraordinarily extended and flat. At > 1 
Gyr after the encounter (bottom right panel of Figure 2) the 
surface density of the ring is ~ 2 orders of magnitude lower 
than in the 'Cartwheel phase', and is almost comparable 
with the density of the disc. The disc now extends up to 
100 — 130 kpc, showing a flat surface density (in a loga- 
rithmic scale). Such flat and huge discs have been observed 
only in GLSBs. 



We compare the properties of our simulated galaxies 
with those of the observed GLSBs. In particular, we con- 
sider a sample of four GLSBs, which have been deeply stud- 
ied: UGC 6614, Malin 1, Malin 2 and NGC 7589. By means 
of the package TIPsMll we derive the surface brightness pro- 
files of the simulated galaxies. We then compare the simu- 
lated surface brightness profiles with the observed surface 
brightness profiles of the four galaxies. Figure 3 shows the 
best-matches between observations and simulations (see Mapelli| 
et al. 2008b for details). The simulations match quite well 
the observations. We also derive the star formation history 
of the simulated galaxies, and we find a good agreement 
with the observed star formation rate, when available. For 
example, the star formation rate inferred from observations 
for Malin 1 is ~ 0.1 yr~^, whereas the simulated star 
formation rate is 0.3 M© yr"^ (Impey & Bothun 1989). 
The simulations also match some interesting morphological 
feature of many observed GLSBs, such as the existence of 
a bar (Pickering et al. 1997; Barth 2007). 

Furthermore, we study the kinematics of the simulated 
galaxies. We rotate the simulated galaxy by the observed in- 
clination angle and then we derive the velocity of gas parti- 
cles along the line-of- sight, in order to produce simulated HI 
spectra. The results (shown in Mapelli et al. 2008b) match 
quite well the data by Pickering et al. (1997). With a similar 
technique we also derive the rotation curves of the simu- 
lated galaxy. In particular, we rotate the simulated galaxy 
by the observed inclination angle, we divide the galaxy into 
concentric annuli and for each of them we calculate the av- 
erage local velocity along the line-of- sight. The results are 
shown in Fig. 4. We make this kind of plots, instead than 
simply show the circular velocity, because we want to do 
something as similar as possible to what the observers do. 
We stress that the rotation curves that we obtain can strongly 
deviate from the circular velocity, as the simulated galaxies 
have strong non-circular motions. In fact, with the method 
that we adopt (and that is adopted in Pickering et al. 1997) 
the entire velocity along the line-of- sight is considered and 
there is no way of distinguishing between circular and ra- 
dial motions. Thus, the rotation curves shown in Fig. 4 gen- 
erally overestimate the velocity with respect to the circular 
velocity. This solves the apparent angular momentum dis- 
crepancy between the circular velocities of ring galaxies and 
of GLSBs. We also note that Pickering et al. (1997) admit 
the existence of strong non-circular motions in the GLSBs 
they analyze, and especially in Malin 1. Thus, it may be im- 
portant to re-analyze the existing data or to take new data 
and to perform a new analysis, in order to measure non- 
circular motions. In our simulations we have information 
about non-circular motions, and we can give predictions that 
may be confirmed (or rejected) by future observations. In 
Fig. 5 the ratio between radial (Vrad) and tangential (vtan) 
velocity is shown, as a function of radius. In all the four sim- 
ulated galaxies there is a clear trend: in the inner region Vrad 
is negligible with respect with vtan, because the central part 
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Fig. 3 i?— band stellar surface brightness profile of the GLSB sample (in units of magnitude per arcsec^). Green crosses: 
data points from Pickering et al. (1997) for UGC 6614, Malin 2 and NGC 7589, and from Moore & Parker (2006) for 
Malin 1. The 1 cr errors are of the same order of magnitude as the points. Red solid line: stellar surface brightness profile in 
R magnitude derived from the simulations. From top to bottom and from left to right: UGC 6614 data and run C of Mapelli 
et al. (2008b) at time t = 0.5 Gyr, Malin 1 data and run C of MapelH et al. (2008b) at t = 1.4 Gyr, Malin 2 data and run B 
of MapelH et al. (2008b) at t = 1.0 Gyr, NGC 7589 and run A of MapelH et al. (2008b) at t = 0.5 Gyr. 



of the galaxy is regularly rotating. In the peripheral regions 
Vrad is similar to vtan, because a part of the ejected matter 
is rapidly falling back to the centre and the remaining is still 
expanding in the outer part of the ring. 



4 Discussion and conclusions 

Our simulations show that collisional ring galaxies in the 
late stages of their evolution (~ 0.5 — 1.5 Gyr) match many 
properties of the GLSBs (surface brightness profile, HI spec-| 
trum, rotation curves, etc.). This is an interesting result, as 
the origin of GLSBs is an open issue, so far. In particular, 
the proposed scenario allows to explain the origin of GLSBs 
within the current cosmological model. Furthermore, the es- 
timated local density of GLSBs (ncLSB) and that of ring 
galaxies (u^q) are comparable: ncLSB ^ 2.7 x 10"'' Mpc~ 
(where h is the Hubble constant in units of 100 km s~^ 
Mpc-^) and urq ^ 5.4 x 10"^ Mpc~^ (Few & Madore 
1986). Ring galaxies are ~ 20 times more numerous than 
GLSBs, a difference which can be explained by the fact that 
magnitude-limited surveys are biased against GLSBs and 
that a fraction of ring galaxies likely does not evolve into 
GLSBs. 

Further theoretical and especially observational tests are 
required, in order to confirm this model of GLSB forma- 
tion. First, kinematic data of GLSBs should be re-analyzed 
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Fig. 4 Rotation curves of the GLSB sample. Green 
crosses are observational data from Pickering et al. (1997). 
1 a errors are shown. Blue open circles connected by the 
solid line are the simulations. The simulations are the same 
as in Fig. 3. 
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Mapelli & Moore: Ring galaxies and GLSBs 
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Fig. 5 Ratio between radial (Vrad) and tangential veloc- 
ity (^tan) of gas in the simulations. The simulations are the 
same as in Fig. 3. 



in order to find possible radial motions. For example, in- 
teresting preliminary results have been found by Coccato 
et al. (2007), who analyze the kinematic data of the GLSB 
ESO 323-G064. Second, it would be interesting to search 
for galaxies which are at the intermediate stage between 
ring galaxies and GLSBs. A possible example is UGC 6614, 
which is considered a GLSB, but shows a ring-like feature 
at ~ 10 — 20 kpc from the centre. In our model this ring- 
like feature may be explained with the secondary ring pro- 
duced by the galaxy interaction. Another interesting object 
is UGC 7069 (Fig. 6), which looks like a typical ring galaxy, 
but has a huge radius (~ 50 — 60 kpc, Ghosh & Mapelli 
2008). These peculiar objects deserve further studies. Fi- 
nally, the redshift evolution of the number of ring galaxies 
and of GLSBs is still an open issue. This can be addressed 
by means of both observations and cosmological simula- 
tions. Interesting predictions for ring galaxies have been re- 
cently presented by D'Onghia, Mapelli & Moore (2008). 
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Fig. 6 Top panel: SDSS composite image of UGC 7069. 
Bottom panel: projected density of stars in a simulation 
which matches the main properties of UGC 7069. See 
Ghosh and MapelH (2008) for details. 



© 2006 WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



www.an-journal.org 



